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rn : ABSTRACT 
O ■ 
O ; 

^ ' We calculate the nuclear composition of matter in accretion disks surrounding 

Q ■ stellar mass black holes as are thought to accompany gamma-ray bursts (GRBs). 

^ ■ We follow a mass element in the accretion disk starting at the point of nuclear 

■ dissociation and calculate the evolution of the electron fraction due to electron, 

^ ■ positron, electron neutrino and electron antineutrino captures. We find that the 

CN ■ neutronization of the disk material by electron capture can be reversed by neu- 

^ ' trino interactions in the inner regions of disks with accretion rates of 1 Mq/s and 

O ■ higher. For these cases the inner disk regions are optically thick to neutrinos, and 

00 ' so to estimate the emitted neutrino fluxes we find the surface of last scattering 

^ ■ for the neutrinos (the equivalent of the proto-neutron star neutrinosphere) for 

O ■ each optically thick disk model. We also estimate the influence of neutrino in- 

^ ■ teractions on the neutron-to-proton ratio in outflows from GRB accretion disks, 
and find it can be significant even when the disk is optically thin to neutrinos. 

o ■ 

^ ■ Subject headings: gamma ray:bursts-nucleosynthesis-accretion disks 

' 1. Introduction 

It is increasingly thought that the progenitor of a gamma-ray burst may be a rapidly 
accreting black hole (see, e.g., Meszaros (2002) for a review), formed from either the collapse 
of a massive star (Woosley 1993; Paczyhski 1998; MacFadyen and Woosley 1999; MacFadyen, 
Woosley, and Heger 2001; Vietri and Stella 1998) or the collision of compact binaries (Eichler 
et al. 1989; Narayan, Paczyhski, and Piran 1992; Ruffert and Janka 1999; Fryer et al. 1999; 
Paczyhski 1991; Ruffert and Janka 2001). Either case is thought to result in the formation 
of a black hole of one to several solar masses surrounded by a debris torus of similar mass 
accreting at rates of fractions to several solar masses per second. Models of GRB accretion 
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disks calculated by MacFadyen and Woosley (1999) and Narayan, Piran, and Kumar (2001) 
suggest that much of this disk material is not actually accreted but is lost in vigorous winds. 
If so, GRB's may be important contributors to local, and possibly galactic, nucleosynthesis, 
especially for rare species such as r-process nuclei. 

The nuclear composition of the ejected material is determined by processing both in 
the disk and in the outflow. Within the disk, the nuclei dissociate as the accreting mate- 
rial becomes hotter and denser closer to the black hole. The nuclear composition is then 
determined by the forward and reverse processes 

e~-\-p n-\-Ve (la) 
-\-n p + Pg. (lb) 

Electron capture on the liberated protons tends to raise the neutron to proton ratio n/p, 
or, equivalently, to reduce the electron fraction Ye ~ 1/(1 + n/p). The extent to which the 
remaining three capture processes mitigate this effect is our primary concern here. 

Accretion disk composition has been previously investigated by Pruet, Woosley, and 
Hoffman (2002). Pruet, Woosley, and Hoffman (2002) dynamically follow the evolution of 
the electron fraction within the disk models of Popham, Woosley, and Fryer (1999) (hereafter 
PWF). However, they consider only the forward reactions of Eqns. la and lb, and so are 
hmited to disk models with lower accretion rates (m < 1, where m = M/{MqS~^)) where 
neutrinos are likely less important. Here, we proceed dynamically as in Pruet, Woosley, and 
Hoffman (2002) but include electron neutrino and antineutrino capture in our evaluation 
of Ye. As anticipated, our solutions for the optically thin PWF disk models do not differ 
significantly from Pruet, Woosley, and Hoffman (2002) with the inclusion of neutrinos. For 
accretion rates of 1 Mq/s and higher, though, neutrinos become quite important in the inner 
regions of the disk. For these higher accretion rates, we shift to the disk models of DiMatteo, 
Perna, and Narayan (2002) (hereafter DPN) that exphcitly include the effects of neutrino 
opacity. In all cases, we also calculate the equilibrium electron fraction for comparison, and 
find it coincides with the dynamical Yg only in the innermost regions of the disk. 

We also model a simple adiabatic outflow to illustrate the importance of neutrino inter- 
actions to the evolution of Ye in this material. We find that while the electron and positron 
capture rates fall off rapidly as the material expands and cools, the neutrino rates remain 
relatively constant. This effect is due entirely to the disk geometry, and so the importance 
of neutrino interactions in the outflow depends quite sensitively on the outflow parame- 
ters. Therefore, the neutrino captures may be important in the outflow, even in the PWF 
scenarios. 

In section 2 we describe the calculation of the electron fraction within the disk. Section 3 
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contains the results of our disk calculations and comparisons with previous work. The simple 
outflow model is discussed in section 4. 



2. Disk Calculation 

We calculate the evolution of Yg in the disk by following a mass element as it moves 
radially from the point at which all nuclei are dissociated (taken here to be T = 10^° K) to 
the inner radius of the disk (just outside the Schwarzschild radius). The method by which our 
calculation proceeds depends on the neutrino opacity of the disk model. The outer regions of 
the disk, regardless of accretion rate, are optically thin to neutrinos - any neutrinos produced 
via electron or positron captm^e typically escape the disk without interacting. Here neutrino 
capture rates are calculated using the neutrino fluxes produced by electron/positron capture, 
and the rates of all four reactions are used to find the change in Yg. In the inner regions of 
the disk, particularly if the accretion rate is 1 Mq/s or higher, the neutrinos become trapped, 
and so the calculation of Yg is not as straightforward. 

The boundary between these zones can be estimated by finding the neutrino optical 
depth as a function of radius within the disk. At a given radius r, we estimate the neutrino 
optical depth as: 

= pK^H = H/l,, (2) 

where p is the baryon density, H is the scale height of the accretion disk, the Ki, is the 
neutrino opacity, and li, is the neutrino mean free path. An equivalent expression holds 
for the antineutrino optical depth. For both neutrinos and antineutrinos, the opacities and 
mean free paths include charged-current and neutral-current neutrino-nucleon interactions 
and neutrino-electron and neutrino- antineutrino scattering. 

As long as Tj,, Tp < 2/3, the disk is said to be optically thin and the calculation proceeds 
as in section 2.1 below. Where t^, Tp > 2/3, both neutrinos and antineutrinos are trapped 
and the capture reactions of Eqns. la and lb come into equilibrium. This is described 
in section 2.3. Between the optically thin and optically thick regions is a zone where the 
neutrinos are 'partially' trapped - > 2/3 while is still < 2/3. This intermediate zone is 
treated in section 2.2. Fig. 1 shows a profile of an accretion disk (a DPN model with m = 10) 
illustrating the three zones defined as above. The region to the right of the long-dashed line 
is optically thin, the region to the left of the short-dashed line is optically thick, and partial 
trapping occurs in between. Using this definition, Eq. 2 for the optical depth, our optical 
depths agree quite well with DPN, differing by a maximum of 20% in the m = 10 case. 

We should note that both sets of disk models we employ (PWF, DPN) assume an 
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electron fraction of 1/2, although in many of the models we find drops far below 1/2. 
While we expect the uncertainty this introduces into our results is small compared to the 
uncertainties in the disk models themselves, both types of physics should eventually be 
included in one self-consistent calculation. 



2.1. Optically Thin Region 

The evolution of in the optically thin regions is calculated by 

dY 

Y.^Vr-^^X^-X^Y, (3) 
dr 

where Vr is the radial velocity of the disk material, Ai = Ai^^n+Ae+n, and A2 = Ai+Ap^p + Ag-p. 
\-p-i ^v^m ^e+m ^-nd Ap^p are the rates for the forward and reverse reactions in Eqns. la and 
lb. Neutron decay is unimportant relative to the short (~ seconds) dynamical timescale and 
so is not included. The above expression for Ye also neglects a (small) general relativistic 
correction and is clearly only valid once all nuclei have dissociated. Except where noted, all 
calculations begin assuming Yp = Y^ = 0.5 at T = 10^'' K. We follow a mass element by 
stepping through the radial zones j, of width rfr^, of the disk model from the point Tq where 
T ~ 10^° K to the inner radius of the disk rj^ax, just outside the Schwarzschild radius. At 
each zone j, we explicitly evolve Ye according to Eqn. 3 above, so that 

Ye,j - Ye,j^i + AYe = FeJ-l + (Al - X2Ye,j-i)drj/Vr,j. (4) 



(5) 



The electron and positron capture rates Xg-p and Ae+„ are given by 

/"OO 

K-p = / (^e{Ee,-Anp)cfe{Ee,IIe)dEe 
poo 

K+n = / Cre{Ee,Ar,p)cfe{Ee,-IJ.e)dEe (6) 

(7) 

where 
and 

Here A„p 1.293 MeV is the neutron-proton mass difference, jie is the electron chemical 
potential, mgC^ ~ 0.511 MeV is the electron mass, gw ~ 0.65 is the dimensionless weak 
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coupling constant, Mwc^ ~ 80.9 GeV is the mass of the W boson, cy ^ I and ca ~ 1-26 
are the vector and axial vector couplings, and 9c ~ 13° is the Cabibbo angle. The electron 
chemical potential is set by the requirement n^- — ne+ = pN^Ye, where ne-,ne+ are the 
Fermi-Dirac number densities for the electrons and positrons, respectively, p is the baryon 
density, and Na is Avogadro's number. 

The neutrino and antineutrino capture rates X^^^n a-nd Ap^p are given by 

POO 

K^n = a / {E + KpfCdE (10) 
Jo 

POO 

Apep = a / {E- Kpf(t)fdE, (11) 

J Anp+meC^ 

where 

and (f)'^, (f)'^ are the effective neutrino and antineutrino fluxes in units of l/(cm^-s-keV). 
Note that here as elsewhere in the paper, we have neglected general and special relativistic 
effects, c.f. Pruet, Fuller, & Cardall (2001). 

The effective fluxes at each radius within the disk include not only the neutrinos and an- 
tineutrinos produced by electron and positron capture at that radius but contributions from 
neutrinos and antineutrinos produced throughout the disk. This introduces a complication 
in our overall disk calculation, as the effective neutrino flux in the outer radial zones depends 
on the flux emitted from the inner zones, which in turn is sensitive to the yet uncalculated 
local Yg. We therefore require several iterations of our overall disk calculation. In the first 
iteration, we dynamically solve for Yg throughout the disk according to Eqn. 4 with the 
neutrino capture rates Xi,^n and Ap^p set to zero. At each zone j, we calculate the neutrino 
and antineutrino fluxes emitted from electron and positron capture within that zone by 

(p^j = ae{Ee,-Anp)cfe{Ee,Pe)npdrj (13) 

(f>9,j = ae{Ee,Anp)cfe{Ee,-pe)nndrj, (14) 

where and (ppj give the number of neutrinos and antineutrinos emitted per second per 
keV per unit area of the emitting region, dvj is the zone thickness, (Xe and arc defined in 
Eqns. 8 and 9 above, and Up and n„ arc the proton and neutron number densities, respectively. 
In the second and subsequent iterations, we repeat the evolution of Y^ through the disk 
including X^^n and Xp^p, with the effective neutrino fluxes at each zone determined by the 
emitted neutrino fluxes of the previous iteration (described below). We continue to iterate 
the disk calculation until Ye throughout the disk differs from the previous iteration by an 
average of less than 1%. 
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The calculation of the effective neutrino fluxes at each zone requires some care, par- 
ticularly since the disk geometry seems to necessitate a fully three-dimensional calculation. 
In general, each disk zone j is a thin cylindrical shell of radius Vj, thickness drj, and total 
height equal to twice the local scale height Hj of the disk. The neutrinos produced in each 
zone are taken to be emitted from the surface area of this shell, and so the emitted fluxes 
(puj and (ppj are in actuality functions of height in the disk, since we expect the temperature 
and density, and therefore the electron and positron capture rates, to drop with height. This 
fact complicates the evaluation of the effective neutrino flux, since (nominally) the effective 
flux in zone i is 

jmax „ 

0f = E / <f>'^,^dn,jA7r, (15) 

3=0 

where d^^ij is the solid angle subtended by zone j as viewed from zone i. Instead of em- 
barking on the involved and numerically expensive exact evaluation of this expression, we 
take the emitted neutrino flux from each zone to be a height-adjusted constant, so that (p^j 
is decoupled from the solid angles dQij. The solid angles, depicted in Fig. 2, can then be 
evaluated separately, as described in the appendix. 

In order to estimate the height adjustment for the emitted neutrino flux, we first model 
the thermodynamics of the disk as a function of height. We simply take the disk to be 
an adiabatic gas in hydrostatic equilibrium, where the total pressure is the sum of gas and 
radiation pressures. We use the temperature and density profiles thus generated to estimate 
the emitted neutrino fiux as a function of height, (j),^j{z), by applying Eqns. 13 and 14 at 
each height z. The correction factor is then 

1-1 • if 4'vj{z)dz 

height correction = hj — , -rj — > (16) 

where is the emitted neutrino flux at z = 0. An equivalent expression holds for the 
antineutrinos. 

The effective fluxes 0^^, (jff require yet another correction. In principle, a neutrino 
diffusion calculation is required here since as we approach the optically thick part of the 
disk, the neutrinos will start to scatter. The transition between the optically thin and thick 
regions is much less sharp than in the protoneutron star of the core collapse supernova. While 
not undertaking the full-blown calculation, we still approximate this scattering by including 
an extinction correction. As the neutrino flux emitted from zone j travels to zone i, it will 
be at least partially extinguished by scattering and captures in the intervening disk material. 
This effect is quite difficult to account for exactly, in particular because the neutrinos can 
follow any number of possible paths from zone j to zone i. We find a minimum extinction 
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correction by assuming the neutrinos follow the minimum path from j to i, so that 

i 

extinction correction = Xij — JJ^ exp(— drfe/Z^^''), (17) 

k=j 

where dvk is the width of zone k and I'j^^ is the neutrino mean free path for captures only 
in zone k. We expect this to be a reasonable approximation as long as > rj. For < rj, 
the maximum possible path length between j and i can greatly exceed the minimum, and 
so Eqn. 17 will be an underestimate. However, since the neutrino flux typically rises sharply 
in the interior of the disk, only a small fraction of the effective neutrino flux in a given zone 
comes from larger radii and so we don't expect the underestimate to be a problem. 

The effective neutrino flux of zone i is therefore 




where hj and Xij are given by Eqns. 16 and 17, respectively, and J dilij/An is calculated as 
described in the appendix. The effective antineutrino flux is calculated in the same fashion. 
If the disk is entirely optically thin {t^, Tp < 2/3 everywhere), then in the second iteration 
and beyond we find the effective neutrino flux as above for each zone, and dynamically 
calculate Yg as previously described. If not, we treat only the optically thin part of the disk 
in this fashion and switch our approach to that described below once trapping sets in. 



2.2. Partial Trapping 

Once Tj, drops below 2/3, the neutrinos begin to be trapped vertically within the disk. 
Again, this would best be described using a one- or two-dimensional Boltzmann neutrino 
transport calculation, particularly since changes rather slowly with radius. However, we 
don't believe this level of sophistication is currently warranted given the still-large uncertain- 
ties in the disk models themselves. Instead, once the neutrinos become trapped as defined 
above, we assume the neutrinos become thermalized and the forward and reverse reactions 
in Eqn. la come into equilibrium. Therefore we replace the effective neutrino flux for these 
zones with a flux calculated from a simple Fermi-Dirac distribution, 

iFD _ 9yC ^ /.„s 

Here g^, is one and /ii, is calculated from the equilibrium condition /i,^ + /i^ — /J'e + l-i'p, or 
III, — lie — A, where ji — Hn — IJ'p is calculated using the Lattimer and Swesty (1991) equation 
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of state. Ti, is taken to be the local temperature of the disk material. It is corrected for 
height in a manner similar to Eqn. 16, except here the integral is over one mean free path 
instead of the entire height of the disk. Aside from this one change, the calculation proceeds 
as in Section 2.1 above. 



2.3. Optically Thick Region 

Once Tp drops below 2/3 as well, both neutrino species are considered trapped and 
we begin solving for Yg assuming the reactions of Eqns. la and lb are in equilibrium. In 
addition, we assume lepton number = + Yy, where Y^ = {n^^ — UpJ/pNA, leaks very 
slowly out of the optically thick region so that Y^ can effectively be considered a constant. 
The value of Yl is set to Ye at the edge of the optically thick region, and we search for values 
of Ye and Yi, at every radial zone j that satisfy 

pLe-+pie+ = (20) 
pLu + pip = (21) 
P'u + P'n = /^e + /^p, (22) 

since the temperature here is well over an MeV. Again, fi = fin — /ip is calculated as in 
Lattimer and Swesty (1991), and the number densities ne-,ne+,n,,, and n^, arc taken to be 
Fermi-Dirac with chemical potentials as above and temperatures equal to the disk tempera- 
ture in that zone. Solving this system of equations proceeds as follows: for each zone j, we 
guess a Yg, find /ig by inverting Ug- — ne+ — pN jCYei calculate fi and solve for //,^, find rii, 
and np, and check to see if the resulting plus Yg keeps Y^ constant. Such a procedure is 
similar to that described in Beloborodov (2003), except that in that calculation there is no 
way to estimate the appropriate value of Yl at the edge of the optically thick region, since 
the system was not considered dynamically. 

In addition to finding Y^, we need to estimate the neutrino flux emitted from the optically 
thick region. In doing so we are guided by the example of the proto- neutron star (PNS). In 
the PNS case, the neutrino luminosities can be found from the rate at which gravitational 
binding energy is released, and the neutrino temperatures can be estimated by determining 
where the neutrinos decouple from the baryons and electrons within the protoneutron star. 
Here we cannot follow this prescription for finding the neutrino luminosities as we have no 
means to determine how much energy is lost into the black hole; instead, we calculate the 
neutrino and antineutrino fluxes from thermal Fermi-Dirac distributions (Eqn. 19) with one 
correction. We add to these fluxes the neutrinos emitted from the optically thin portions 
of the disk into the optically thick region. This small correction is necessary to ensure that 
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these neutrinos don't simply disappear into the optically thick region. This has a negligible 

effect on the result and if we had a perfect blackbody, it wouldn't be needed. Contrary to 
the case of the protoneutron star, the boundary between trapping and no-trapping is not 
at all sharp. We in effect artificially draw such a boundary, and therefore use this and the 
extinction correction to try to correct for this sharp boundary approximation. We should 
note that the above prescription for the neutrino and antineutrino fluxes gives values of the 
same order of magnitude as those calculated from energy considerations assuming no energy 
is lost into the black hole. 

In order to follow the PNS example to evaluate the neutrino temperatures, we need to 
translate the concept of a neutrinosphere to the disk geometry. In general, the last scattering 
surface for a sphere can be defined as the radius Ri, at which the neutrino optical depth Ti, 
becomes 2/3: 



where is the neutrino opacity and p is the baryon density. It makes little sense to apply 
this expression to the plane of the disk (integrating from the outer edge of the disk inward), 
as the neutrinos that, by the above definition, arc 'trapped' horizontally can easily escape 
the disk vertically. Instead, we adopt Eqn. 23 to calculate where the neutrinos decouple 
vertically as a function of disk radius. We take the vertical disk structure determined as 
described in section 2.1, and calculate the height at which the following expression is 
satisfied: 



where hmax » H and the neutrino opacities Ky are taken to be equal to {lyp) ^, where ly 
is the neutrino mean free path, here averaged over energy. 

The shapes of the neutrino decoupling surfaces thus calculated are shown in Figs. 3 
and 4. The long dashed lines show the decoupling heights hy as a function of radius for the 
neutrinos; the short dashed lines are for the antineutrinos. The scale heights H (solid line) are 
plotted for comparison. Contrary to the PNS case, these are not spherical. As Fig. 3 shows, 
the decoupling heights hy in the DPN rh — 1 model steadily increase with decreasing radius, 
resulting in wedge-shaped decouphng surfaces. In the DPN m — 10 model, the optically 
thick region is so large that the decoupling surfaces are additionally shaped by the physical 
height of the disk. Within the radial extent of the optically thick region the scale height H 
decreases appreciably. Therefore while hy/H and the emerging neutrino and antineutrino 
temperatures continue to increase with decreasing radius, the decoupling heights hy actually 
decrease, resulting in the rounded 'neutrinosurf aces' seen in Fig. 4. 




(23) 




(24) 



A further simplification is made to facilitate the calculation of the neutrino and an- 
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tineutrino fluxes emitted horizontally into the outer disk. Instead of the complicated shapes 
shown in Figs. 3 and 4, we take the 'neutrinosurface' to be a cylinder, with radius equal 
to the outer tip of the 'neutrinosurface' (where hy 0) and height equal to the maximum 
decoupling height. The temperature and chemical potential of the emerging neutrinos are 
taken to be their local values at r^. The temperatures that we find range from around 
Ty = 2.5 — 4.5 MeV and Tp = 3.6 — 5.1 MeV depending on the model. The neutrino and 
antineutrino fluxes found in this way are included in the effective neutrino flux calculations 
described in section 2.1. 



3. Results - Ye in Disk 

The results from our calculations of the electron fraction in the disk models of PWF 
and DPN are summarized in Table 1, which contains the final values of Y^ in the innermost 
radial zones of each disk model. In this table rh is the mass accretion rate, a is the black 
hole spin parameter, and a is the viscosity. 

For disks with accretion rates m < 0.1, the evolution of the electron fraction is domi- 
nated by electron capture. For these low accretion rates the entire disk is typically optically 
thin to neutrinos, and so neutrino interactions play a limited role in the disk. For disks 
with higher accretion rates, m > 1, Y^ is set by all four capture reactions. In these disks, 
electron capture initially drives Ye to very low values (Ye <0.1). As the mass element spirals 
inward, the material becomes optically thick to neutrinos and neutrino interactions become 
increasingly important. Neutrino and positron capture significantly raise Ye in the inner disk 
to the final values in Table 1. 



3.1. m<l 

Figs. 5 and 6 show two representative calculations of the evolution of Ye in optically thin 
accretion disks, using PWF disk models with m = 0.1, alpha viscosity a = 0.1, and black 
hole spin parameter a = (Fig. 5) and a = 0.95 (Fig. 6). For each model, the solid line is 
our full calculation and the dotted line is our calculation without the neutrino interactions. 
The dashed lines show for comparison the equilibrium electron fractions 

= ^, (25) 

where Ai and A2 are defined as in Eqn. 3. 

The calculated electron fraction shown in Fig. 5 is particularly representative of the 
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PWF m — 0.01, 0.1 models. The temperature and density are relatively low so we see only 
some electron capture and a small drop in Y^, halted by positron capture within r ~ 100 km. 
The capture rates are too slow for Yg to equilibrate, and neutrinos have almost no effect, as 
stated in Pruet, Woosley, and Hoffman (2002). Models with lower alpha viscosity a < 0.1 
have higher densities and so more electron capture (and a correspondingly lower Ye), but 
neutrinos are similarly unimportant. However, in the high spin (a > 0) models the neutrino 
capture reactions begin to play a role. In these models the portion of the disk close to the 
Schwarzschild radius is hotter and many times denser than the equivalent disk with a — 0. 
In these conditions the neutrinos may even become trapped. 

An example of this is shown in Fig. 6. As in Fig. 5, Y^ drops initially due to electron 
capture, and then rises slightly as positron capture becomes more important. The steep 
drop in Ye beginning at r ~ 30 km is due to a rapid increase in density (raising the electron 
capture rate) combined with a small dip in temperature (dropping the positron density and 
therefore the positron capture rate). Once the temperature again begins to rise, positron 
capture rebounds and Ye levels off. Within a small portion (10 < r < 14 km) of this 
very dense and hot region the neutrinos become trapped. The corresponding jump in the 
neutrino capture rate raises Yg within this narrow zone. In this small region the crossing 
time is comparable to neutrino capture time, about half a millisecond. However, this is 
not a concern for our calculation, since we do not assume weak equilibrium. The fact that 
the neutrinos (not the antineutrinos) are trapped means that their flux and spectra are 
calculated assuming thermal equilibrium (see section 2.2), but Ye is calculated in this region 
by integration of the weak rates. 

Neutrinos have comparatively little effect elsewhere within the disk, as can be seen 
by comparing the full calculation in Fig. 6 (solid line) to the calculation without neutrino 
interactions (dotted line). It is not immediately apparent that this should be the 
it may seem that the high flux from the optically thick region should produce appreciable 
neutrino capture throughout the disk. This does not happen because the effective neutrino 
flux from the optically thick region (or, in fact, from any region in the disk) is rapidly 
extinguished via two separate effects. The first is the extinction correction from section 2.1. 
The neutrino opacities in the regions just outside of the optically thick zone are still quite 
high, so a portion of the flux is lost to neutrino capture. The second, and most important, 
effect is entirely geometric. At any radius, much of the neutrino and antineutrino flux will 
leave the top or bottom of the disk without interacting with disk material. As a result, the 
effective flux drops off much faster than, for example, in the proto-neutron star case, where 
the flux falls as r~^. This is illustrated in Fig. 7, which shows the effective neutrino flux from 
the optically thick region over the net flux emitted from that region as a function of r (solid 
line). The dashed line is the same quantity with the extinction correction removed, so that 
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only the geometric effect is reducing tfie flux. The dotted line shows r ^ for comparison. 

In addition, in this model the neutrino and antineutrino fluxes effectively negate each 
other's influence in the outer disk. Even though the neutrino flux from the center of the disk 
is much higher than the antineutrino flux, the neutrino opacities are everywhere larger than 
for the antineutrinos, leading to greater extinction of the neutrino flux as it passes through 
the disk. This effect is also shown in Fig. 7, where the effective flux over the emitted flux 
for antineutrinos is given by the dot-dashed line. As a result of this, the effective neutrino 
flux for r > 25 km is only slightly higher than the effective antineutrino flux, and the net 
impact of the neutrinos is small. The extent that this effect influences Ye is illustrated by 
the dot-dashed line in Fig. 6. It shows a calculation where the antineutrino interactions only 
are removed and the neutrinos, acting alone, more appreciably raise for r > 25 km. Still, 
the geometric effect dominates the extinction of the neutrino flux and so even here Y^ is not 
radically altered outside of the optically thick region. 



3.2. m>l 

Both neutrinos and antineutrinos become trapped in the inner regions of accretion disks 
with higher accretion rates (m > 1). For our calculations of Y^. in such disks we would prefer 
to use disk models that incorporate neutrino transport effects. We therefore switch to the 
DPN disk models, which, unlike the PWF models, self-consistently include the effects of 
neutrino opacity. We find that given the PWF parameters, their model is optically thin, 
whereas given the DPN parameters their model is optically thick. Our results for Y^ in the 
inner regions of these disks are very sensitive to this choice. 

Fig. 8 shows our calculated electron fraction for the m — 1 DPN model. Again, the 
solid line is our full calculation, the dotted line is our calculation with the neutrinos removed, 
and the dashed line is the equilibrium electron fraction. Also included for comparison is our 
full calculation for the PWF m = 1 model with a = 0.1 and a = 0. The temperatures and 
densities here are significantly higher than for m = 0.1, and, as a result, electron capture 
quickly drives Y^ to a very low value. Positron capture moderates this drop within r ~ 250 
km. As shown, neutrino interactions radically alter Y^ in the DPN model. Once neutrino 
trapping sets in at r ~ 70 km, thermalization decreases the neutrino temperature slightly 
but raises the effective flux; the latter wins out and the neutrino capture rates increase 
markedly. This drives Ye back up to 0.26 before antineutrino trapping sets in and relowers 
Ye to ~ 0.24. 

As in the PWF m = 0.1, a — 0.95 disk model from section 3.1, we see that neutrinos 
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have a noticeable impact within the optically thick region and much less influence outside 
of that region. Here the sections of the disk where neutrinos arc at least partially trapped 
are larger and so the resulting in the inner disk (0.25 compared to closer to 0.1 in the 
PWF model or the DPN model without neutrinos) is more dramatically altered. in the 
partially trapped region is particularly sensitive to the neutrino physics; when we artificially 
modify T^, by only a few percent we find the resulting Y^ in this region can change by 30% 
or more. However, outside of the trapping regions, the influence of neutrinos is limited by 
the same factors: extinction, competition between neutrino and antineutrino captures, and, 
most importantly, the geometric effect. 

Fig. 9a shows the net neutrino (solid lines) and antineutrino (dashed lines) fluxes as a 
function of radius for the DPN rh — 1 model. The thin lines show the net fluxes emitted in 
the optically thin regions, and the flUed and unflUed circles show the net neutrino and an- 
tineutrino fluxes, respectively, emitted from the optically thick regions. The thick lines show 
the calculated effective fluxes at each radius due to contributions from the rest of the disk 
(except for the 'partially trapped' region, 38 < r < 68 km, where the effective neutrino flux 
is from a Fermi-Dirac distribution as described in section 2.2). It is important to note that in 
the optically thin regions the total effective fluxes at each radius are significantly higher than 
the emitted fluxes. In these regions the effective fluxes are dominated by contributions from 
the rest of the disk, particularly from the optically thin regions interior to that zone where 
the emitted fluxes are higher and from the optically thick region. Additionally, though the 
emitted antineutrino flux is everywhere smaller than the neutrino flux (and often orders of 
magnitude smaller), the effective antineutrino flux is less than a factor of two smaller than 
the effective neutrino flux for most of the disk with r > 100 km. 

Fig. 10 shows our calculated electron fraction for the m = 10 DPN model. As in Fig. 8, 
the solid line is our full calculation, the dotted line is our calculation with the neutrinos 
removed (for the optically thin and partial trapping zones only), and the dashed line is the 
equilibrium Ye. The evolution of Y^ here is similar to the m = 1 case, with one notable 
exception. The inner regions of the disk are much hotter than the m = 1 model, and 
so positron capture plays a much larger role. This is indicated in the calculation without 
neutrinos, where Yg begins to increase at r ~ 300 km. The greatest influence, though, is in 
the optically thick region (r < 160 km). Here Y^ continues to increase, albeit at a slower 
rate, even when the antineutrinos are trapped. Since Yl is (approximately) constant in 
this region, becomes negative, which favors antineutrinos over neutrinos. As a result the 
antineutrino flux from the optically thick region is signiflcantly higher and more energetic 
than the neutrino flux. 

The emitted and effective neutrino fluxes in this model are shown in Fig. 9b, where the 
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lines and points are defined as in Fig. 9a. Again, the major difference between the two panels 
is the higher antineutrino flux from the optically thick region. This carries over to the outer 
regions of the disk, where the effective antineutrino flux is approximately three times that 
of the neutrino flux. Still, this difference isn't large, and the effective flux suffers from the 
same geometric effect as the other models. As a result, the influence of the neutrinos within 
the disk is again largely confined to the regions in and around where they are trapped. 

The results of this section are summarized in Figs. 11 and 12. These figures show the 
electron fraction Y^. as a function of position for inspiraling mass elements in the DPN disk 
models fn = 1 and 10, respectively. We note the importance of positron, neutrino, and 
antineutrino capture by marking the points where each interaction, when omitted, changes 
Ye by more than 10%. 

We have calculated the fluxes coming everywhere from the surface of the disk and they 
affect many important pieces of physics. They alter the electron fraction in material which 
flows off the disk and they also determine the neutrino-antineutrino annhilation rate above 
the disk. A preliminary analysis of the former is in the next section. An analysis of the 
latter will be done in future work, but to give a point of reference: for the DPN model with 
m — 1, along the z-axis, at a height of 8 x lO^cm (again neglecting relativistic corrections) 
the neutrino number flux is 1.27 x lO'^^cm"^, the antineutrino number flux is 5.0 x lO'^^cm^^. 
The average energies of the neutrinos and antineutrinos at this point arc 15.1 MeV and 15.7 
MeV respectively, where as the rms energies are 16.4 MeV and 17.1 MeV respectively. 



We further examine the influence of neutrinos on nucleosynthesis in GRB's by following 
material from the disk as it is ejected in a wind. Neutrinos leaving the disk will interact with 
the outflow material, thus changing its electron fraction. Our goal here is not to develop a 
realistic outflow model but to estimate the import of neutrino interactions on nucleosynthesis 
in the wind. To this end we choose a simple constant velocity, adiabatic outflow. We assume 
a constant mass outflow rate proportional to 47rr^pt', which gives p ~ r"^, where r here is 
the outflow radius in the spherical approximation. The temperature in the ejecta is found 
from the density and the entropy per baryon. Here we assume a constant entropy per baryon 
in the outflow. The entropy per baryon includes the contributions from relativistic particles 
and nucleons, as in Qian and Woosley (1996): 



4. Preliminary Outflow Model 




(26) 
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where TmbV is the temperature in MeV and pio is the density in units of lO^'^ g/cm^. We 
start with a mass element in the disk, launch it with velocity v equal to the escape velocity, 
and follow the evolution of the electron fraction in the ejecta as described in section 2.1 until 
the temperature drops to 10^° K. 

An important difference between this calculation and that of section 2.1 is the evaluation 

of the effective neutrino and antincutrino fluxes at each point. We simplify the evaluation of 
the effective fluxes for points above the disk by assuming that the disk is essentially flat as 
viewed from above. Therefore each disk zone is no longer a cylindrical shell but a flat ring 
with radius rj and width drj. The emitted fluxes from each region first need to be adjusted for 
this change in geometry Instead of multiplying the number of neutrinos emitted per second 
per keV per volume {a^cfe from Eqn. 13) by the width dvj of the zone for a cylindrical 
emitting surface, we multiply this quantity by the depth dj of the emitting zone. For the 
totally optically thin regions, the depth is just the total height of the disk, dj — 2H. For 
zones where the neutrino mean free path l„ is less than 2H, we have dj = ly. The height 
correction hj discussed in section 2.1 is also applicable here. The emitted neutrino flux 
from each zone j is therefore 



An equivalent expression holds for the antineutrinos. 

The emitted fluxes from the optically thick region are found using the 'neutrinosurfaccs' 
calculated as described in section 2.3. At each radius within the optically thick region, the 
emitted fluxes are calculated from Fermi-Dirac distributions according to Eqn. 19. The 
temperature in this expression is taken to be equal to the vertical disk temperature at the 
appropriate decoupling height hy or hy for that radius, and the chemical potential is equal 
to iiy or lip at that disk radius. These temperatures vary along the 'neutrinosurface'. For 



example for DPN m = 10, 2.4 < < 6.0 and 3.6 < < 6.8, while for DPN m = 1, 
4.5 <Ty< 5.3 and 5.1 < Tp < 5.4. 



The effective fluxes are found by integrating the emitted fluxes over the entire disk, 
similar to in section 2.1 but with very different geometry. Here the appropriate solid angle 
is not so easily decoupled from the emitted fluxes, so we evaluate the full integral for the 
effective flux ^'J^ at each point {x,y,z) in the outflow. 




(27) 
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In the above expression, (f)'^^{Ey,9,(f)) is the emitted flux 0'^^ at Vj — [2;^tan^^ + + 
2x2; tan ^ cos 0]^/^, and ^max(0) is the maximum value for 9 at the disk edge. This angu- 
lar geometry is illustrated in Fig. 13. 
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Fig. 14 shows the minimum variation in Ye for two such calculations for outflows from 
the DPN m — 10 disk model. The solid line is for an outflow from close to the black 
hole, at Tdisk ~ 35 km. The short-dashed line is the equilibrium Y^, for the same outflow. 
Electron capture and antineutrino capture here combine to lower Kg from its initial value 
in the disk. As the electron capture rate drops, the equilibrium Ye (meaning that Ye that 
would obtain if the weak rates were in equilibrium) is increasingly set by neutrino and 
antineutrino capture. However, Y^. itself levels off before then since the very rapid outflow 
velocity {v — Vesc — 1-5 x 10^° cm/s) causes Ye to quickly fall out of equilibrium. The second 
outflow example, the long-dashed hne in Fig. 14, starts from the disk region just outside of 
the optically thick zone and immediately falls out of equilibrium. Here the reason is not so 
much that the outflow velocity is high, but that the equilibrium Ye rises sharply above the 
disk. This is entirely a neutrino-driven effect. Within the disk, the neutrino and antineutrino 
fluxes are quickly diminished by extinction and, most importantly, the geometric effect, so 
the equilibrium Ye is set largely by electron and positron capture. Above the disk, however, 
the outflow material is exposed to the high neutrino and antineutrino fluxes leaving the disk, 
and so the jump in the neutrino and antineutrino capture rates readjusts the equilibrium 
Ye, here to a much higher value. In this example, the temperature falls below 10^° K before 
Ye levels off, but still the outflow velocity is too great (here v — Vesc — 5.7 x 10^ cm/s) for 
equilibrium to be established. 

Figs. 15 and 16 further illustrate the importance of neutrino capture in the wind. Fig. 15 
compares the outflow Yg calculated as above to that calculated without including neutrino 
and antineutrino captures. The percent change in Ye, (Y^^^" - Y^"") x 100%/^^"°'', is 
plotted for outflows from every ten radial zones in the DPN disk model with m — 10. In 
outflows from the outer disk, the neutrino flux is small and so is the percent change in Ye. 
Above the inner disk, particularly above the optically thick region, the neutrino flux is much 
higher and Ye is influenced accordingly. The reason the neutrinos can have such a large 
influence is shown in Fig. 16. It plots the capture rates in two outflow models, (a) DPN 
m— 10 and (b) PWF m = 0.1, a — 0.95, as a function of height. In both cases, the electron 
and positron capture rates fall quickly, dropping below the neutrino and antineutrino rates 
within a few hundred kilometers. The neutrino rates drop much more slowly due to the 
geometry of the disk. If the disk emitted neutrinos uniformly and the outflow remained 
above the disk, the rates would be almost constant. In our physical disks, this doesn't hold 
exactly since most of the neutrinos are from the inner regions of the disk, but still the rates 
fall much more slowly than, for example, ~ as in the proto-neutron star case. 

It is important to note that our outflow model only conservatively estimates the influence 
of neutrino interactions in the wind. For example, if the velocity of the outflow is at any 
point slower than the escape velocity at the starting point of the flow, the wind's exposure 
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to the neutrino flux will be correspondingly greater. In addition, we only follow the outflow 
until the temperature drops below 10^° K; we expect neutrino interactions to continue to 
play an important role in the subsequent nucleosynthesis. This is true even for disks that 
are almost entirely optically thin (see, for example. Fig. 16b). 

5. Conclusions 

Here we have made a detailed analysis of the importance of the weak rates in the disks 

of gamma ray bursts. We follow mass elements as they spiral into the center of these disks, 
keeping track of all the rates, the emitted neutrino and antineutrino fluxes, and the electron 
fraction. We have included not only the effect of electron and positron capture but also of 
neutrino and antineutrino capture. Neutrino and antineutrino capture play an essential role 
for models with high accretion rates, where the neutrinos become trapped. Here the electron 
fraction can change by factors of several over calculations where these rates are not included. 
In addition, neutrinos also play an important role in models with lower accretion rates if 
they have high black hole spin parameter. 

Even if the weak rates come into equilibrium toward the center of the disk, it is still 
necessary to follow the complete evolution of the mass element as we have done. This is 
because the lepton number of the material as it falls into the optically thick regions goes 
into determining the electron fraction and fluxes of neutrinos in these areas. 

As part of this study, we have calculated the 'neutrinosurf aces' or surfaces that define 
the regions that are optically think to neutrinos. We find that these are not spheres as in 
the protoneutron star case, but take on different shapes, usually more wedgelike. 

As part of the larger picture, we wish to determine the nucleosynthesis that will come 
from the mass outflow from these disks. This will involve a hydrodynamic calculation, in 
addition to simply knowing the Y^. on the disk which we have calculated here. However, 
whatever this outflow looks like, it is certain that the neutrino and antineutrino flux from 
the disk will hit it from behind as it moves out. This will change the electron fraction and 
therefore also any calculation of nucleosynthesis. We have made an estimate of the minimum 
influence of this effect by using a constant velocity for the outffowing material which is equal 
to the escape velocity. Even in this case, the neutrinos change the Yg in the outflow by 
10%-60%. 

Future studies of nucleosynthesis from these disks will produce exciting results. What- 
ever the results will be, it is certain that the neutrino interactions will have to be included 
in the calculations. 
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A. Appendix 

Here we describe the evaluation of the integral over the solid angle J dQij from Eqn. 18, 
where dQi^j is the solid angle subtended by zone j as viewed from zone i within the disk. 
As depicted in Fig. 2, we take zone j to be a uniform cylindrical shell with radius rj, height 
2Hj, and thickness drj and the viewpoint in zone i to be a point at the midplane of the disk 
a distance n from the black hole. The general form of the integral is as follows: 



dnij= / / sm0d0d(f). (Al) 

J<Pa Jea{4>) 



The solid angle geometry is dependent upon whether zone j is interior or exterior to zone i, 
and so the limits of integration for each case are treated separately below. Once the limits 
are determined we solve the integrals in Eqn. Al numerically using the extended Simpson's 
rule. 

The case where zone j is exterior to zone i is pictured in more detail in Fig. 17. For 
convenience we arrange the coordinate clXGS clS shown, with the origin at the viewpoint in 
zone 2, the -\-x axis running along the midplane of the disk and through the black hole, and 
the z axis extending vertically above and below the disk. In every case where Vj > r^, part 
of zone j is blocked by the event horizon and possibly an optically thick region, as depicted 
in the top view detail of Fig. 17. Therefore the limits of the integral over the coordinate 
become 

0a = sin-(^^) (A2) 

06 = 27r - sin- (^^^) , (A3) 

where Rblock is the Schwarzschild radius or, if present, the radius of the optically thick 
region. 

The limits of the integral over 9, as illustrated in the side view detail of Fig. 17, are 

= tan-^(^^^ (A4) 
^6(0) = 7r-^„(0), (A5) 

where r(0) is the distance from the origin to the midplane of zone j as a function of (p. The 
distance r(0) is found from the law of cosines 

— r1 + — 2rir cos (A6) 
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to be 

r(0) = n COS + (r| - sin^ 0) . (A7) 

The case where zone j is interior to zone i is illustrated in Fig. 18. The limits of the (J) 
integral, as shown in the top view detail, are 

,-1 / ^'J 



06 = +sin-^(^-^J. (A9) 

The side view detail of Fig. 18 shows the limits of the integral over 6. The expressions 
for 9a{(f>) and 9b{(f>) are identical to Eqns. A4 and A5, except here the distance r(0) is the 
negative solution of Eqn. A6, 

r (0) = Ti cos - (rj - r'^ sin^ 0)^/^ (AlO) 
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Table 1. The Electron Fraction near the Event Horizon 
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Fig. 1. — Disk profile for tlie DPN m = 10 disk model. The scale height H of the disk is 
plotted versus disk radius r, from the Schwarzschild radius to where the temperature drops 
below 10^*^ K. The vertical long-dashed hue marks the radius where the neutrino optical 
depth Ti, drops below 2/3, and the short-dashed line marks the equivalent point for the 
antineutrinos. 
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Fig. 2. — Illustration of the solid angle, d^lij, subtended by zone j as viewed from zone i, 
for where j is interior to i, i.e. rj < rj, and where j is exterior to i, i.e. rj > Vi. Note that for 
Tj > Ti, part of zone j is blocked by the black hole and (possibly) the optically thick region. 
We calculate the flux at the point n from every accessible place on the surface at Vj. 
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Fig. 3. — Decoupling heights hi, for neutrinos (long dashes) and antineutrinos (short dashes) 
versus radius in the DPN disk model m = 1.0. The sohd hne marks the scale height H of 
the disk. 
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Fig. 4. — Same as Fig. 3 for the DPN disk model m — 10. 
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Fig. 5. — Electron fraction Kg as a function of radius for the PWF models with m = 0.1, 
alpha viscosity a = 0.1, and black hole spin parameter a = 0. The solid line shows from 
our full calculation, while the dotted line is our calculation without neutrino interactions. 
The dashed line is the that would obtain, were the system in weak equilibrium. 




Fig. 6. — Identical to Fig. 5 for the PWF disk model with m = 0.1, a = 0.1, and a = 0.95. 
Again the solid line is the full calculation and the dotted is the calculation without neutrinos. 
The rise at the center of the disk is due to neutrino trapping. The additional dot-dashed 
line shows our calculation with just the antineutrinos turned off. 
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Fig. 7. — Effective neutrino fluxes from the optically thick region divided by the net flux 
emitted by the optically thick region, for the PWF disk model with m = 0.1, a = 0.1, and 
a — 0.95. The dot-dashed hne is for the antineutrino flux and the solid and dotted hues 
are for the neutrino flux calculated with and without the extinction correction, respectively. 
The dashed line shows for comparison. 
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Fig. 8. — Electron fraction as a function of radius for the DPN and PWF disk models 
with rh — 1, black hole spin parameter a — 0, and alpha viscosity a — 0.1. The dot-dashed 
hne shows Ye calculated with the PWF disk model, while the remaining three lines are from 
the DPN disk model. The sohd line shows Yg from our full calculation, the dotted hne is for 
a calculation without neutrino interactions, and the dashed line shows the equilibrium Yg. 
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Fig. 9. — Net neutrino fluxes versus radius for the DPN disk models with (a) m = 1 and 
(b) rh — 10. The dark sohd (dashed) hues show the net effective (anti) neutrino flux at each 
radius, while the lighter weight lines show the net neutrino flux emitted from each optically 
thin zone. The fllled and unfllled circles depict the net antineutrino and neutrino fluxes, 
respectively, from the optically thick region and are plotted at the effective 'neutrinosurface' 
radii r^,- 
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Fig. 10. — Electron fraction as a, function of radius for the DPN disk model with m = 10, 
a — 0, and a = 0.1. As in Fig. 8, the solid line is the full calculation with neutrino 
interactions, the dotted line without, and the dashed line shows the equilibrium Ye. 
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Fig. 11. — Electron fraction Yf, as a function of disk position {x,y) for an inspiraling mass 
element in the DPN disk model with m = 1.0, a = 0, and a = 0.1. The markings indicate 
where each capture reaction becomes important, defined as where the absence of that reaction 
in the calculation affects Yg by greater than 10%. 
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Fig. 12. — Same as Fig. 11 for DPN disk model m = 10. 
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Fig. 13. — Illustration of the angular geometry relevant to calculating the effective neutrino 
flux for a point {x, y, z) above the disk. 
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Fig. 14. — Electron fraction versus height in the outflow originating from two locations 
within the DPN disk model m = 10: close to the black hole (r ~ 35 km, solid line) and at 
the outer edge of the optically thick region (r ~ 250 km, long-dashed line). The short-dashed 
line and dotted line show the equilibrium Yg for each outflow. 
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Fig. 15.— Percent change in Y^, (y/'^^" - Y^'"") x 100/^^'"°'^, is plotted for outflows from 
DPN disk model m = 10. Y^^^^'^ includes the effects of neutrino interactions, while Y^"" is 
calculated with the neutrino and antineutrino capture rates set to zero. 
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Fig. 16. — Capture rates versus height in the outflow originating from two disk models, (a) 
DPN rh — 10 and (b) PWF m = 0.1, a — 0.95. For both models, the solid line is the electron 
capture rate, the long-dashed line is the positron capture rate, the short-dashed line is the 
antineutrino capture rate, and the dotted line is the neutrino capture rate. As shown, the 
neutrino and antineutrino capture rates remain high well above the disk. 
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Fig. 17. — More detailed illustration of the solid angle, dQ.i^j^ subtended by zone j where j 
is exterior to i, i.e. Vj > Vi. 



-39- 




side view 




Fig. 18. — More detailed illustration of the solid angle, d^ij, subtended by zone j as viewed 
from zone i, for where j is interior to i, i.e. rj < ri. 



